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Abstract. The static n-body problem of General Relativity states 
that there are, under a reasonable energy condition, no static n- 
body configurations for n > 1, provided the configuration of the 
bodies satisfies a suitable separation condition. In this paper we 
solve this problem in the case that there exists a closed, noncom- 
pact, totally geodesic surface disjoint from the bodies. This covers 
the situation where the configuration has a refiection symmetry 
across a noncompact surface disjoint from the bodies. 



1. Introduction and background 

A classical result in Newtonian gravity is that there can be no static 
n-body configuration for which the bodies are separated by a plane dis- 
joint from the bodies. On the other hand one can concoct static 2-body 
configurations in Newtonian theory [BS] with both bodies being con- 
tractible and one body sufficiently non-convex so that the convex hulls 
of the bodies intersect. Analogous configurations exist for relativistic 
bodies (work in progress by L. Andersson, the first author, and B. G. 
Schmidt). For n > 1 and assuming a suitable energy condition, it is 
reasonable to conjecture a relativistic analogue of the Newtonian result 
stated above; that is, n-body static configurations should be impossi- 
ble provided some separation condition for the bodies is satisfied. The 
work [Mu] has some results on the static n-body conjecture, but no 
theorem under easily stated conditions. In the present paper we show 
(see Theorem 2.2) that an asymptotically fiat triple (M, V, g) with non- 
negative scalar curvature which is static vacuum outside a compact set 
and in a neighborhood of a closed, embedded, noncompact, totally ge- 
odesic surface is trivial. This solves the static n-body problem in the 
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case that the configuration has a refiection symmetry across a noncom- 
pact surface which is disjoint from the matter regions (see Theorem 
2.3). 

Recall that static spacctimes are 4-manifolds with a metric of Lorentz 
signature which have a Killing vector field which is complete, every- 
where timelike, and hypersurface orthogonal. General Relativity stud- 
ies such spacetimes subject to the Einstein equations G^v = 8nGTf^,j 
(see [W]). Such solutions describe the gravitational fields of time in- 
dependent, non-rotating sources. Static spacetimes can be written as 
warped products R x M with metric ds'^ of the form 

ds'^ = -V'^{x) df + gij{x) dx'dx^ (1.1) 

with V a positive function and g a Riemannian metric on the 3-manifold 
M. The Einstein equations then take the form 

AV = A7rGV{p + r) (1.2) 

and 

VRij - DiDjV = AnGV [{p - r) g^j + 2Tij] , (1.3) 

where p and Tij = T(^ij) are respectively the energy density and the stress 
tensor in the rest system of the matter and r = Ti is the trace. We are 
interested in solutions to these equations corresponding to n isolated 
bodies. By this we mean the following: First the 3-manifold (M, g) is 
asymptotically flat with V tending to 1 at infinity. Secondly we assume 
that the support of the matter fields p, Tij is contained in n disjoint 
compact connected sets Vt^, with Vt^ open with smooth boundary dVl^ 
for r = 1,... ,n. Finally we assume that all fields are sufficiently 
smooth (even analytic) except across dVl^r where p, and the normal 
components of d'^gtj, d^V will in general have jump discontinuities. 
We require also that g and V be G^ across the boundaries. Let us 
remark that taking the trace of (1.3) and using (1.2) we recover the 
time symmetric initial value constraint 

R = IGnGp (1.4) 

and taking a divergence of (1.3), using (1.2) and the contracted Bianchi 
identity, we find that 

Dj(Vr/)+pDiV^O, (1.5) 

which plays the role of equilibrium condition for the matter variables. 
In order for this condition to hold distributionally across the boundaries 
we require the additional boundary condition 

Ti^rij\gn^ = (1.6) 
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that is, the stress should have zero normal components to the boundary 
of the bodies. In many models of continuum mechanics the stress 
tensor is a functional of a collection of matter fields and their first 
derivatives, which renders equation (1.5) a quasilinear second order 
PDE with Neumann-type boundary conditions (1.6). For perfect fiuids 
one has that Tij = pgij with p > in fir and p a given positive non- 
decreasing function of p in M"*". There are different energy conditions 
which one might impose on the matter variables (see [HE]), the weakest 
one being that p > 0, which is sufficient for the positive mass theorem 
[SY] to be vahd. Finally one might mention here the case of black holes, 
in which the regions Ur^lr are missing, but instead at the boundaries 
VldUr = with dflr being totally geodesic surfaces. 

Historically, the 'no-body situation', i.e. n = 0, implies that (M, V, g) 
is trivial (Minkowski) in the sense that V = 1 and (M, g) is fiat was 
the first to be classified. This is the content of a classical result in [L] 
if M is assumed to be diffeomorphic to (the proof extends easily 
to all topologies). After many partial results it was recently shown 
by Masood-ul-Alam [Ma] that when matter is composed of a perfect 
fiuid we must have n — 1 and the spacetime is spherically symmetric; in 
particular, Schwarzschild in the vacuum region. These spherical models 
have been studied extensively [HRU]. Solutions for n = 1 without 
(spatial) symmetries, for sources composed of ideally elastic material 
have been constructed in [ABS]. For black holes it is known that n has 
to be 1 and the solution is isometric to the exterior of a Schwarzschild 
black hole. This has been shown in [BM] in the nondegenerate case 
and in [C] generally. 

2. Separating surfaces 

Let (M, g) be an asymptotically fiat Riemannian three manifold. 
We allow the possibility that M has a finite mimbcr q > 1 of ends 
Ma, 1 < a < Q', each being asymptotically fiat. Recall that the static 
vacuum equations are given by VRij — Vij — and AV = for a 
positive function V where Rij denotes the Ricci tensor of g and Vij the 
covariant hessian of V taken with respect to g. We will be interested 
here in metrics which are static vacuum solutions outside a compact 
set, and at the very least have nonnegative scalar curvature everywhere. 

We will consider a surface S which is noncompact, connected and 
properly embedded in M. We first show that if such a surface is to- 
tally geodesic, then it has a finite number of ends each of which is 
asymptotic to a plane in one of the asymptotically fiat ends of M at 
infinity. Precisely we show that there is a compact subset K oi M such 
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that for each a, Ma H {S \ K) is equal to a finite union of graphs of 
functions fp, 1 < p < ka, over a Euchdean plane (in suitable coordi- 
nates on Mq,) such that fp approaches a constant and its derivatives 
decay at an appropriate rate. 

Proposition 2.1. Let S be a noncompact, connected, totally geodesic 
surface properly embedded in M. Assume that S fl is unbounded 
in the end Ma- There exist asymptotically fiat coordinates defined on 
Ma such that outside a compact set K the surface S fl {Ma \ K) is the 
union of ka > i graphs of functions = fp{x^,x'^) for 1 < p < ka 
such that there are constants ap so that fp — ap decays like 1/r' and 
the derivatives of the fp decay correspondingly faster, where r' = r'^ = 



^J{x^Y + [x'^y. Note that this description holds for each of the ends 
Ma for which S fl Ma is unbounded and the number ka depends on a 
as do the coordinates and the functions fp. (We take ka — if S H Ma 
is bounded.) We omit the dependence of the coordinates and the fp on 
a for notational convenience. 

Moreover, for a sufficiently large the compact subset of S given by 
Sa — S n {K U (U^^ilr^ < (t})) is a compact surface with boundary 
curve Co- (having k — X]q=i ka components) such that the Euler char- 
acteristic x{.Sa) is equal to x{S) and linio-^oo Jc^ k ds = 2T:k where k 
is the geodesic curvature of the oriented curve C„ in S. 

Proof. Our argument will work separately on each end, so throughout 
we focus attention on one end Ma such that S fl Ma is unbounded and 
we omit explicit reference to a unless needed for clarity Prom the work 
of [B2] there exist coordinates on Ma defined outside a compact set K 
such that g is equal to a Schwarzschild metric up to order r~^, that is 



where m is the ADM mass. (We use the notation 0{r~^) to denote a 
term which is bounded by a constant times r^^ and whose derivatives 
up to a fixed order decay correspondingly faster.) Since S is embedded 
and the manifold Ma \ K may be chosen to be simply connected (for 
example we can take it to be diffeomorphic to the exterior of a ball 
in M?) it follows that S is orientable. We choose the orientation for 
M and hence for S determined by the coordinates Ob ^ Ob ^ Ob , and let 
Ci and 62 be an oriented local orthonormal basis for S relative to the 
metric g. It then follows that the length N of the 2- vector ei A 62 with 
respect to the Euclidean metric is 1 + 0(r^^). Therefore using the fact 
that S is totally geodesic with respect to g we have D^Ji^ei A 62)] = 
for q; = 1, 2. Letting V denote the Euclidean connection, observe that 
the difference tensor T — D — V is ol order since it is given in 




^,, = (l + 2m/r)5,, + 0(r-^) 
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Euclidean coordinates by the Christoffel symbols of g, so we have 

= Ve„(ei Aesj+TeJeiAes). 

Prom this we see that Ve„(ei A 62) is 0(r~^) and therefore 

Ve^N = 7V-^(Ve„(ei A 62)) • (ei A 62) = 0{r-^). 

Now the length of the second fundamental form of S with respect to 
the Euclidean metric is the Euclidean magnitude of V(A^~^ei A 62) 
taken along S (since A^~^ei A 62 is the Euclidean unit tangent plane), 
and therefore the length of the Euchdean second fundamental form is 
0(r-2). 

Note: The argument above shows that if g = 6 + 0(r~^), then the 
magnitudes of the second fundamental form of S taken with respect to 
the indicated metrics satisfy the inequality \As\ < c\Ag\ +cr^^ since in 
this case the difference tensor is 0(r~^). 

Let (To be a radius to be chosen large, and let M„ o- denote the part of 
Ma exterior to the open ball of radius a > aQ. Let Sq > and consider 
the rescaled surface S{ao) = eo/ao{S n M„,<,) C \ B^^{0). The 
length of the second fundamental form of S'(cro) is then equal to ctq/so 
times that of S at corresponding points, and distances arc changed by a 
factor of Eq/ctq, so we see that the second fundamental form of S'(o"o) at 
a point X is bounded by c(£o/c"o)|2^|~^- Since S is connected, we see that 
either S{ao) has a single component without boundary or it has ka> 1 
components 5'p((To), 1 < p < ka, each with boundary on ^^^^(O). In 
the former case it follows from Proposition 3.1 (next section) that for 
ctq sufficiently large (hence the second fundamental form small with 
quadratic decay), S is the graph of a function / over a plane which 
we may take to be the x^x^-plane, and that the second derivatives of 
/ decay like 0((r')^^), and the first derivatives like 0((r')~^). In the 
second case Proposition 3.1 implies that each of the Sp{aQ) may be so 
described as the graph of a function fp with the same decay conditions. 
Note that since S is embedded each of the Sp{ao) is a graph over the 
same plane. 

Scaling back to the original surface S we obtain the description of 
Sr\{Ma\K) as a union of graphs. To get the required decay, we use the 
Schwarzschild form of the 1/r term in the metric expansion. We observe 
that the metric g defined by ^ = (1 + m/r)'''^g has the property that 
g — 5 + 0{r~'^). Using the well known relation for second fundamental 
forms of conformally related metrics we see 



= + (1 + m/r)-^v{l + m/r)g 
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where z/ denotes the unit normal of S with respect to g and for a 
function (p, we use i>{(p) to denote the derivative of (p in the direction 
{>. Since Ag = and from the asymptotic behavior of the fp we see that 
on the graph of fp we have z/ is plus or minus + 0(r~^), so we have 
\Ag\ = {\/3m\x^\/r^) + 0(r^'^). From the fact that first derivatives 
of / decay like 0((r')~^) it follows that fp is bounded by O(logr'). 
Putting — fp in the bound on the second fundamental form, we 
see that \A^\ = O ( (log r)r~^). Since the metric g is Euclidean up to 
terms of order r~^, we use the Note above to improve the decay on the 
Euclidean second fundamental form to O ( (log r)r~^). This can then be 
used to show that fp is bounded and has a limit Up at infinity. Putting 
this information back into the second fundamental form bound tells us 
finally that the second derivatives of fp decay like 0((r')~^), and this 
implies the desired asymptotic decay. 

The final statement on the behavior of the total geodesic curvature 
follows from the easily checked fact that the geodesic curvature of Co- 
is equal to l/cr + 0{a~^) while the length of each component of Ca is 
equal to 27rcT + 0(1). □ 

Theorem 2.2. Assume that M is static vacuum outside a compact 

set and has R> everywhere. Suppose there is a closed, noncompact, 
totally geodesic surface S such thai g is static vacuum, in a neighborhood 
of S . It follows that M is isometric to the Euclidean space M^. 

Proof. Let V be the static potential defined in a neighborhood of S and 
outside a compact set of M. We first show that V is identically 1 on 
S and that S is flat (zero Gauss curvature). To see this, we choose a 
local orthonormal frame so that the are tangential for a = 1, 2 and 
63 is normal to S. We then take the tangential trace of (1.3) to obtain 

VRaa = Ka = AgV 

where we have used that fact that S is totally geodesic to write the 
trace of the covariant derivatives on M in terms of the intrinsic Laplace 
operator on S. (It would be sufficient here that S be minimal.) Now 
the Gauss equation tells us that since S is totally geodesic we have 

Raa Rafiafi 

where K is the intrinsic Gauss curvature of the surface S. Since R = 
in the static vacuum region due to (1.2), this implies that R33 = —Raa, 
and therefore Raa — K. Thus we see that the restriction of V to 
S satisfies the equation /S.sV — KV — 0. Now we let be as in 
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Proposition 2.1, and apply the Gauss-Bonnet theorem to obtain 
/ K da^ 27rx(-S) - / k ds. 



The totally geodesic condition implies that K — -R1212 is bounded by 
a constant times r~^, and thus by Proposition 2.1, K is an integrable 
function on S. Thus we may let a tend to infinity to conclude Jg K da = 
2tix{S) — 2nk < since k > 1 and the Euler characteristic of any 
connected noncompact surface is at most 1. On the other hand we 
have K — V~^AsV, so we may also write 



where u is the outer unit normal along Cp. Since V tends to 1 and 
the derivatives of V decay at least as fast as it follows that the 
boundary term goes to as p goes to infinity and we have 



We therefore conclude that the integral on the right is and hence V 
is constant on S. It follows that V = 1 on S, and from the equation 
satisfied for V that K = on S. It follows moreover that x(<S') = 1, 
and hence S is isometric to the Euclidean M^. 

Now it is a known asymptotic property of the static equations ( [Bl] , [B2] ) , 
that there is a constant m so that 



and that m is equal to the ADM mass. Thus we have shown that m 
is zero, so it follows from the Positive Mass Theorem [SY] that M is 
isometric to the Euchdean M^. This completes the proof. □ 

The following result is a consequence of Theorem 2.2. 

Theorem 2.3. A nontrivial relativistic static n-body configuration can- 
not have a reflection symmetry across a noncompact surface which is 
disjoint from the bodies. 

Proof. Assume we had such a configuration with S being the surface 
fixed by the symmetry F. It would then follow that S is totally geodesic 
since a geodesic a beginning at a point of S and initially tangent to 
S must remain in S since F o o" is a geodesic with the same initial 
conditions and is therefore identical to a. The result now follows from 
Theorem 2.2. 







□ 
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3. A technical result for surfaces in 

In this section we prove the technical result used in the proof of 
Proposition 2.1. That result is the following. 

Proposition 3.1. Assume that S is a dosed, connected, noncompact, 
embedded surface in \ B^^ where Br denotes the closed ball of radius 
r centered at the origin. Assume also that for any point x & S we 
have \A\{x) < cSqIxI'"^ where A denotes the second fundamental form 
of S. If Eq and 6q are sufficiently small, then there exist Euclidean 
coordinates x^, x^, x^ so that any connected component of S (1 (M^ \ Bi) 
is contained in the graph of a function x^ — f{x^,x^) defined for r' — 

{x^Y + {pc^y > 1/2 such that the first and second derivatives of f 
satisfy \df \ < c{r')-^ and \d'^f \ < c(r')-^. 

Proof. We first consider the case in which S fl dB^Q = 0. In this case, 
S is a closed embedded surface in R^. Let P G S* be a point nearest 
the origin and note that \P\ > Eq. We choose Euclidean coordinates 
y^-iU^-iy^ so that P is at the origin and so that ^{P) = where u 
denotes the unit normal vector field to S. There is a neighborhood of 
in iS which is the graph of a function — fi{y^, y^) defined for p' — 

^Jiv^Y + {y'^Y <Rso that < 1. We show that the set of R with 
this property consists of all positive real numbers, and thus the entire 
surface S may be so represented. To see this, let R be the largest radius 
for which such a representation is possible, and use the fundamental 
theorem of calculus along the ray ^{t) — {ty^, ty'^, fi{ty^, ^y^)) to write 

Since \dfi\ < 1 it follows that \^'{t)\ < V^p', and thus we have 

Hy\y\My\y')) - ^1 ^ ^p' £ my\ty'j,{ty\ty'))\ dt. 

Now Ity^jty"^, fi{ty^,ty'^))\ > tp', and thus from the second fundamen- 
tal form bound we have li^iy^jy"^, fi{y^,y'^)) — ^\ < c5o{p')^^. It 
follows that if 60 is chosen sufficiently small we have \df{y^,y'^)\ < 1/2 
for p' < R. This contradicts the choice of R as the largest radius for 
which \df\ < 1. This shows that S is globally given as the graph of 
a function with gradient bounded by 1. Therefore from the second 
fundamental form bound we have |(9^/i| < c5o(p')^^- follows by in- 
tegration as above that the first partials of /i converge to constants 
at infinity, and thus we may change coordinates to x^,x'^, x^ so that S 
is given as x^ — f{x^,x'^) and so that the first derivatives decay hke 
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(r') ^. This gives tlie desired conclusion under the assumption that 
S n dB,, = 0. 

Let us now assume that S fl dB^a 7^ 0- We first analyze the points 
of S which lie on the unit sphere. Let P G S* fl dBi and suppose that 
the tangent plane of 5* at P does not intersect P2£o- If ^0 is sufficiently 
small this implies that a large neighborhood of P on lies arbitrarily 
close to the tangent plane, and hence does not intersect B^^. In this 
case the argument above implies that a connected component of S is 
a global graph and hence we must have been in the first case. There- 
fore it follows that the tangent plane to 5" at P intersects B2eo, 
therefore since £0 is arbitrarily small, i^(P) is arbitrarily close to being 
tangent to the unit sphere. It follows from this that S intersects dBi 
transversally, and that the curves of intersection have small geodesic 
curvature. Since the curve of intersection is embedded, we can see by 
elementary geometry that it must consist of a finite number of curves 
all of which lie in a small neighborhood of a great circle with each curve 
being close to the great circle. 

Now if we consider a point P on one of these curves 7, then we 
choose coordinates y^,y^,y^ so that the point P is (1,0,0) and that 
i/(P) = A neighborhood of P in 5" may then be represented by the 

graph = fi{y^,y'^) with /i of small norm defined over a disk of 
radius 7/8 centered at (1,0). This representation then extends to cover 
a neighborhood of the curve 7 by the graph y^ = fi{y^, y"^) defined for 
1/4 < p' < 3/2. If we now consider the largest value of R for which 
this representation extends to the set 1/4 < p' < P with \dfi\ < 1, 
then we may repeat the argument above to show that R = 00, and 
thus each of the intersection curves lies on a connected component of 
5'n(M^\Pi) which has the required description as a graph of a function 
over the region r' > 1/2 in the plane. Note that the 1/4 is replaced 
by 1/2 since we need to do a slight rotation of coordinates to make 
the tangent plane at infinity to be the x^x'^ -plane. We could replace 
1/2 by any fixed small radius Vq by taking Sq and Sq sufficiently small. 
Since S is embedded, these planes must be parallel, so the description 
holds simultaneously for all components in a fixed system of Euclidean 
coordinates. This completes the proof. 

□ 

References 

[ABS] Andersson, L., Beig, R., and Schmidt, B. G.: Static self-gravitating elastic 
bodies in Einstein gravity, Comm. Pure Appl. Math. 61, 988-1023 (2008). 

[Bl] Beig, R.: Arnowitt-Deser-Misner energy and goo, Phys. Lett. A 69, 153-155 
(1978/79). 



10 



ROBERT BEIG AND RICHARD M. SCHOEN 



[B2] Bcig, R.: The static gravitational field near spatial infinity, Gen Relativity 

Gravitation 12, 439-451 (1980). 
[BS] Beig, R. and Schmidt, B. G.: Celestial mechanics of elastic bodies, Math. 

Z. 258, 381-394 (2008). 
[BM] Bunting, G. L. and Masood-ul-Alam, A. K. M.: Nonexistence of multiple 

black holes in asymptotically Euclidean static vacuum space-time, Gen. 

Relativity Gravitation 19, 147-154 (1987). 
[C] Chrusciel, P.T.: The classification of static vacuum space-times containing 

an asymptotically flat spacclike hypersurface with compact interior, Class. 

Qu. Grav. 16 661-687 (1999). 
[HE] Hawking, S. W. and Ellis, G. F. R., The large scale structure of space- 
time, Cambridge Monographs on Mathematical Physics, No. 1, Cambridge 

University Press, London, 1973. 
[HRU] Hcinzlc, J. M., Rohr, N., and Uggla, C: Dynamical systems approach 

to relativistic spherically symmetric static perfect fluid models. Classical 

Quantum Gravity 20, 4567-4586 (2003). 
[L] Lichnerowicz, A.: Theories relativistes de la gravitation et de 

I'electromagnetisme. relativite generale et theories unitaires, Masson et Cie, 

Paris, 1955. 

[Ma] Masood-ul-Alam, A. K. M.: Proof that stellar models are spherical, Gen. 

Relativity Gravitation 39, 55-85 (2007). 
[Mu] Miiller zum Hagcn, H.: The static two body problem, Proc. Cambridge 

Philos. Soc. 75, 249-260 (1974). 
[SY] Schoen, R. and Yau, S. T.: On the proof of the positive mass conjecture in 

general relativity. Comm. Math. Phys. 65 45-76 (1979). 
[W] Wald, R.: General relativity, University of Chicago Press, Chicago, IL, 1984. 

R. Beig, Faculty of Physics, Gravity Group, University of Vienna, 
Vienna 

R. Schoen, Department of Mathematics, Stanford University, Stan- 
ford, CA 94305 



